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We study the epidemic spreading process following contact dynamics with heavy-tailed waiting 
time distributions. We show both analytically and numerically that the temporal heterogeneity of 
contact dynamics can significantly suppress the disease's transmissibility, hence the size of epidemic 
outbreak, obstructing the spreading process. Furthermore, when the temporal heterogeneity is 
strong enough, one obtains the vanishing transmissibility, hence the lack of epidemic outbreaks for 
any finite recovery time, the condition of which was derived. 



Introduction — Throughout history, epidemics have 
had major influence not only on an individual's health 
but also on global human history [T], proving itself of 
great intellectual and practical interest [2HZ]- The his- 
tory of epidemics [T] shows that most of infectious dis- 
eases leading to pandemic state such as pest, cholera, 
malaria, and SARS (severe acute respiratory syndrome) 
were transmitted by broadcasting through carriers like 
animals and insects or media like polluted water or air 
without requiring direct physical contacts between sus- 
ceptible and infected individuals. On the contrary, dis- 
eases requiring direct contacts for transmission such as 
sexually-transmitted diseases [8j have rarely been able to 
achieve pandemic states, with only a few contemporary 
examples such as AIDS. 

Differences in transmission mechanisms for the two 
classes of infectious diseases may leave strong imprints 
on the characteristics of corresponding disease progres- 
sion dynamics. For diseases transmitted through broad- 
cast, explicit contact dynamics between individuals may 
become less important, and the classical theoretical ap- 
proach based on the assumption of Poisson process and 
well-mixed population could reasonably be applicable [21- 
[7]. Diseases requiring direct contacts for transmission, 
however, are directly governed by patterns of human ac- 
tivity, for a wide range of which from email exchanges 
to sexual contacts are known to exhibit heavy-tailed dy- 
namics [3 HO] , rendering the Poisson approach question- 
able [11]. 

Timings of individual's actions are typically highly het- 
erogeneous. Human activity is dominated by a few bursts 
of activities, with extended periods of quiescence in be- 
tween [3]. Such temporal heterogeneity [T^] can be dic- 
tated by the waiting time distribution P{t), where the 
waiting time t is the time interval between two consec- 
utive actions, and presents a new layer of complexity in 
social dynamics, parallel with the structural heterogene- 
ity or network heterogeneity dictated by the complex net- 
work structure of contacts [13] . Indeed, it has been shown 
that such temporal heterogeneity can significantly affect 
to slow down spreading processes in networks [TT ][T4Hl8] . 
In this paper, we investigate further the effect of tempo- 
ral heterogeneity of contact dynamics on the large-scale 



properties of epidemics. 

The main result of this paper is to show analytically 
that the temporal heterogeneity can significantly impede 
the spreading of epidemics, in stark contrast with the 
network heterogeneity that facilitates epidemic spread- 
ing [3]. Most dramatically, the epidemic outbreak can 
even be completely suppressed, given sufficiently strong 
temporal heterogeneity. We demonstrate this analyti- 
cally by applying renewal theory 'W,'^ to a prototypical 
epidemic model, the susceptible-infected-recovered (SIR) 
model. We derive expressions for the transmissibility and 
thus epidemic threshold for heavy-tailed P(t), to show 
that the epidemic threshold increases with heterogeneity 
of contact dynamics without bound. The analytical pre- 
dictions are well supported by extensive numerical simu- 
lations on random and scale-free networks. We conclude 
the paper by examining the effect of finite timescale and 
with discussion on the implications of our results. 

SIR model with arbitrary contact dynamics — The SIR- 
type model we consider in the paper is formulated as 
follows. Each node is in one of three states, suscepti- 
ble, infected, or recovered. The disease is transmitted 
through the contact between an infected node and its 
susceptible neighbor. In classical approach [5], the con- 
tact dynamics is assumed to be a Poisson process. Here 
we relax such Poisson assumption, and consider that the 
contact dynamics follows independent renewal process 
with a general inter-event time (waiting time) distribu- 
tion P{t). Along the way, the infected node can recover 
autonomously, after when it does not participate in epi- 
demic dynamics. Classically the recovery dynamics is 
also assumed to be Poissonian. In this work we consider 
a fixed recovery time A, to focus on the effect of heavy- 
tailed contact dynamics. 

Key quantities for disease spreading dynamics are the 
so-called transmissibility T and the secondary reproduc- 
tive number i? [US]. T is the probability that an infected 
individual would transmit disease to a susceptible neigh- 
bor before it recovers, and R is the expected number of 
secondary infections per each infected node. Given P{t), 
X controls the "infectiousness" of the disease (and thus 
T). If A is large (small), there is more (less) chance for 
secondary infections. As the whole population becomes 
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FIG. 1. Transmissibility T for (a) the power-law and (b) the 
lognormal P(r) with A = 10"^ (sohd), lO" (dashed), and 10^ 
(dotted), plotted as a function of the power-law exponent a 
and the standard deviations a of P{t), respectively. Tmin = 1 
in (a). The dash-dotted line in (b) denotes the asymptotic 
formula, Eq. (5), for the lognormal Pij). (Insets) Epidemic 
threshold Ac with k = 10, vs. a (a) and a (b). 



composed of either susceptible or recovered nodes in the 
stationary state {t — > oo), the average fraction of recov- 
ered nodes pao in the stationary state measures the "size" 
of epidemic outbreak. One defines the epidemic thresh- 
old Ac for the epidemic outbreak to be the infimum of A 
such that Pao > 0. 

Transmissibility and epidemic threshold — Following 
renewal theory [20] . the transmissibility for contact dy- 
namics following a renewal process with independent, 
identically distributed P{t) and fixed recovery time A 
can be obtained as 

/•oo /"OO /"OO 

T = / 5(A) / S{tR - X)dtRdA = 1 - / g{A)dA. 

(1) 

Here 5(A) is the so-called generation time distribu- 
tion |19j . the distribution of time intervals between the 
moment of infection and the first following contact ac- 
tivity, which in this case is given by the residual waiting 
time distribution, g{A) — P{T)dT, where (t) is 



the mean waiting time. The integral with respect to tji 
accounts for the probability that the node does not re- 
cover during the interval A. The cases with general re- 
covery time distribution p(i^) can be treated by replacing 
the delta function with p(tji) in Eq. (1). 

In Fig. 1, we show the transmissibility T calculated 
for two heavy-tailed distributions that are widely used 
to model bursty dynamics [9l [21] , (i) the power-law dis- 
tribution with exponent a and minimum waiting time 
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for T > To and P{t) = otherwise, and (ii) the lognormal 
distribution with unit mean and variance cr^. 
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In both cases, the transmissibility T decreases with the 
contact dynamics' heterogeneity, dictated by either the 
power-law exponent a or the variance of P{t), and it 
even vanishes as a approaches to 2 (Fig. la) or as the 
variance diverges (Fig. lb). This result clearly demon- 
strates that the heavy-tailed contact dynamics can dras- 
tically suppress, and even completely block, the epidemic 
spreading. 

For a power-law P{t), Eq. (2), the transmissibility can 
be calculated explicitly. It reads 



(a-2) A 



TpL{X;a) = 
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A < r„ 
A > r„ 



(4) 



Evidently, Tp^ decreases with as a decreases and van- 
ishes as TpL ^ (a— 2) as a — )■ 2. For general P{t), T may 
not always be obtained in a simple form. Yet its asymp- 
totic behavior can be more accessible for many cases. 
For example, for the log-normal waiting time distribu- 
tion, Eq. (3), the residual waiting time distribution is 

- In A-i ln(l+CT^) 



obtained as gL^iA; a) 



1 + erf 



where erf(x) denotes the error function. Using properties 
of the error function, one obtains the leading asymptotic 
behavior of T for large a as 
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(5) 



vanishing algebraically with a (Fig. lb). 

Once T is obtained, the epidemic threshold can be 
readily obtained for the process on uncorrelated, tree- 
like networks, by mapping to a branching process |22j . 
From the criticality condition of the branching process, 
the condition for the epidemic outbreak is written as 
R = Tk > 1, where k is the average branching num- 
ber, given by the expected number of neighbors of an 



3 



infected node excluding the parent node. For an uncor- 
related network k is given by the expected remaining 
degree of a node reached by following a randomly chosen 
link, that is K = Y.k^k- l)kP{k) / {k) = {{P) - {k))/{k). 
For power-law P(r), Eq. (2), the epidemic threshold Ac 
is explicitly obtained as 

''''^ 1 r^in [{a - 1) (1 - i)]"'/^""'^ (A, > w). 



(6) 

Ac increases and diverges as a — ?► 2 (Fig. la, inset). 
Therefore, for a = 2 the epidemic outbreak cannot take 
place for any finite A, implying that only unrecoverable 
diseases (A = 00) can spread through the population. 
Similarly, Ac diverges with a for the lognormal -P(t) 
(Fig. lb, inset). Finally, to obtain p^o, one can apply the 
generating-function method [1] based on the mapping to 
bond percolation in which each bond is randomly occu- 
pied with probability T. 

Numerical simulations — We test the validity of the 
analytical predictions on two random network models, 
the Erdos-Renyi (ER) random graphs and the static 
model of scale- free (SF) graphs The numerical sim- 
ulation runs as follows. Initially all nodes are suscepti- 
ble except for one infected node, chosen at random, as 
a seed node. Each connected pair of nodes makes con- 
tacts following a renewal process with the waiting time 
distribution P(t), except for the first contact, the tim- 
ing of which follows the residual waiting time distribu- 
tion. Whenever an infected agent make a contact with 
a susceptible neighbor, the disease spreads, turning the 
susceptible node into infected. Along the way, each in- 
fected node recover after a fixed recovery period, A. The 
process proceeds until there remains no infected agents 
in the network, and the final fraction of recovered nodes 
S is measured. The ensemble-averaged value of 5* over 
independent runs gives the expected outbreak size poo- 

On ER networks, we show the numerical simulation 
results with both the power-law and lognormal P(t), to- 
gether with the theoretical curves (Fig. 2a, b). The the- 
oretical predictions are in excellent agreement with the 
numerical simulations. The epidemic outbreak size con- 
sistently decreases and the epidemic threshold diverges 
with the strength of temporal heterogeneity of contact 
dynamics, dictated by a approaching 2 (Fig. 2a) or di- 
verging a (Fig. 2b). 

On SF networks with asymptotic power-law degree dis- 
tribution Pd{k) ^ k^"' , it is well-known that the epidemic 
spreading is facilitated to the extent that the epidemic 
threshold vanishes in the limit of infinite network size 
when 7 < 3, as k diverges with TV [3]. To verify the im- 
pact of temporal heterogeneity in such a case, we perform 
numerical simulations with power-law P(t) on the SF 
network with 7 = 2.5 (Fig. 2c). The epidemic outbreak 
size decreases as a decreases, meaning that the tempo- 
ral heterogeneity still hinders epidemic spreading in SF 
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FIG. 2. Plots of the final density of recovered nodes (the 
expected outbreak size) poo vs. the recovery time A, of the 
SIR models, (a, b) Results on ER networks with mean degree 
(k) — 10 and A'' = 10^, for (a) the power-law and (b) the 
lognormal -P(t). (c) Results on SF networks with 7 = 2.5 and 
TV = 10" with the power-law -P(t). Symbols denote numerical 
simulation results and lines in (a, b) denotes the theoretical 
curves, in excellent agreement with each other. (Inset) Inverse 
of the numerically simulated poo with A = 10 (rescaled by 
its value p^for Ns = 10'^), plotted against the number of 
samplings of generation times TVs. For a — 2.1 it remains 
constant, whereas it increases logarithmically for a = 2.0 and 
1.9 (same symbols as in main panel). Therefore, as TVs ^ 00, 
Poo is expected to vanish for a < 2, as predicted by the theory. 



networks. As long as a > 2, however, Ac « 0, that is, 
the epidemic outbreak occurs for any nonzero A. In this 
sense, the network heterogeneity dominates over the tem- 
poral heterogeneity, when a > 2. For a < 2, however, 
the temporal heterogeneity can dominate over network 
heterogeneity, suppressing the epidemic outbreak com- 
pletely even for SF networks. In numerical simulation. 
Poo is obtained to be nonzero, albeit small, for a = 2 
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FIG. 3. The epidemic threshold Ac vs. the maximum waiting 
time Tnmx in the power-law waiting time distribution with 
various exponent a, obtained from numerical simulations of 
the SIR model on the scale-free network with 7 = 2.5 and 
= 10^. Dotted lines have slope 1.0 (top) and 0.2 (bottom), 
shown for comparison with the theoretical prediction, Eq. (8). 
Two regimes are separated by the minimum r, Tmin = 1, 
indicated by the shade. 

(and even for a = 1.9), which is to be attributed to the 
finite number of samplings for the time to first contact 
from .g(A) when a < 2. Indeed, poo is found to decay as 
the number of samplings is increased, and thus expected 
to vanish in the infinite-time limit even for SF networks, 
as predicted by the theory (Fig. 2c, inset). 

Effect of the finite cutoff timescale — So far, we have 
assumed that there is no cutoff in the maximum wait- 
ing time in P{t). In reality, however, contact dynamics 
mediating the spreading process takes place over a finite 
time window, bounded, for example, by an individual's 
lifespan. We examine the effect of such a cutoff timescale 
set by the maximum waiting time on the epidemic out- 
break. With the cutoff waiting time Tmax, the generation 
time distribution is given by g{A) = /X"" P{T)dT. 
Let us take P{t) to be a power law with exponent a 
in the range (rmin, Tmax)- For a > 2, Tmax plays only a 
minor effect in the transmissibility, neglectable for large 
T'max- Whereas for 1 < a < 2, reportedly corresponding 
to a number of human activities [TD] , Tmax becomes the 
dominant term in T(A, a) as 

T(\-r„'\ / •^/'''max "* (A < Tmin), /y\ 

^ ' ' I (A/w)^-" (A>r„un). ^ ' 

Finally, the epidemic threshold Ac depends on Tmax as 

> J '^max ^c,P (A < Tmin), 

Ac ~ <^ l/(2-a) Q > ^ . X (8) 

I 'max^cP ^ 'inm^- 

where Ac,p denotes the epidemic threshold for Poisson 
contact dynamics (exponential P{t)). The predicted de- 
pendence of Ac on Tniax is Well supported by the nu- 
merical simulations (Fig. 3). This result shows that the 



more heavy-tailed (smaller a) the contact dynamics is, 
the larger is the impact of long but finite waiting times. 
In the r,nax — > 00 limit, Ac diverges with Tmax, again 
completely suppressing the epidemic outbreak. 

Summary — To summarize, we have shown both ana- 
lytically and numerically that epidemic outbreaks of the 
SIR model can be strongly suppressed, and even com- 
pletely blocked, by the heavy-tailed contact dynamics. 
Applying renewal theory, we have derived the transmissi- 
bility T and epidemic threshold Ac for contact dynamics 
following power-law and lognormal waiting time distri- 
butions. It is shown explicitly that T vanishes (conse- 
quently, Ac diverges) as a — > 2 or cr — > 00, respectively, 
which are specific instances of the general condition for 
diverging Ac, given by g{A)dA — 1 for any finite A. 
As such, temporal heterogeneity is found to exert op- 
posite effect to network heterogeneity, thus competing 
each other. Finally, it is noteworthy that although we 
have specifically formulated our analysis with the SIR 
dynamics with fixed recovery time, the main result of 
suppressing effect of temporal heterogeneity would ap- 
ply under more general epidemic scenarios. For exam- 
ple, we observed similar effect for the SIR dynamics with 
distributed recovery times [^S] and the susceptible- 
infected-susceptible model [3l El [T^ . 

As various kinds of human activities, including sex- 
ual activities [27j, exhibit strongly heavy-tailed dynam- 
ics, our results should bear broad implications on epi- 
demiological understanding and control of, for example, 
sexually transmitted diseases. Another interesting hy- 
pothesis that can be derived from our results is that the 
infectious diseases requiring direct physical contacts for 
transmission should necessarily have extremely long re- 
covery time, were it to accomplish an epidemic outbreak. 
This may offer a Darwinian selection-based explanation 
on why such diseases as AIDS tend to be so resilient to 
cure that it takes excessively long for complete recovery, 
often persisting even throughout an individual's lifetime. 
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